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New precision Cosmic Microwave Background Radiation (CMBR) anisotropy data are beginning 
to constrain physics beyond the standard model, for example in the form of additional light particle 
species. These constraints are complementary to what can be obtained from big bang nucleosynthesis 
(BBN) considerations because they apply to much later times. We derive a constraint on the 
equivalent number of neutrino species, N v , from the presently available data. Specifically we analyse 
two different CMBR data sets to test the robustness of our results. Analyzing only CMBR data 
yields an upper bound of N„ < 17 (95% confidence). Adding large scale structure (LSS) data from 
the PSC-z survey tightens the upper bound slightly. However, the addition of LSS data gives a 
non-trivial lower bound of N v > 1.5/2.5 (95% confidence) for the two data sets. This is the first 
independent indication of the presence of the cosmological neutrino background which is predicted 
by the standard model, and seen in big bang nucleosynthesis. The value N v — is disfavoured at 
3a and 4a for the two data sets respectively. 

PACS numbers: 98.70.Vc, 14.60. St, 13.35. Hb 



I. INTRODUCTION 

Precision measurements of the anisotropy in the Cos- 
mic Microwave Background Radiation (CMBR) have re- 
cently begun to probe cosmology with high precision. 
The measurements have delivered remarkably strong sup- 
port for inflation for the standard inflationary paradigm, 
i.e. a flat geometry and a initial fluctuation power spec- 
trum which is close to scale invariant. Because of the 
high precision of the current measurement it is also pos- 
sible to probe various other parameters of the standard 
model. In the present paper we study the current limits 
on the relativistic energy density during recombination. 
The energy density is usually parameterized in terms of 
N v , the equivalent number of standard model neutrino 
species 



(1) 



The standard model prediction is N v ~ 3.04, where the 
0.04 comes from the fact that neutrinos are not com- 
pletely decoupled during the electron-positron annihila- 
tion in the early universe |Q . 

BBN considerations give the bound M 



2 < A^.bbn < 4 (95% confidence) 
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A bound on this parameter has been derived previously 
from CMBR data However, it was pointed out 

by Kneller et al. |7j] that the bound is quite sensitive to 
assumptions about other cosmological parameters. 

In the present paper we discuss in detail degeneracies 
between N v and various other cosmological parameters, 
particularly the Hubble parameter, Hq. Using two dif- 
ferent compiled data sets we derive bounds on N v . We 



then go on to discuss the influence of including data from 
large scale structure surveys. It turns out that including 
LSS data significantly narrows the allowed region for N v . 

Apart from providing a fairly robust upper limit on N Vl 
the main result of the present paper is that the cosmic 
neutrino background has been detected at more than the 
3a level (i.e. N v = is disallowed at the 99.7% level). 
The presence of the neutrino background is also detected 
by big bang nucleosynthesis data (see Eq. (||)). However, 
this is the first independent cosmological detection. The 
standard model value N v = 3 is in all cases within 2a of 
the maximum of the likelihood function, so there is no 
evidence for deviations from the standard model in the 
present data. 



II. CMBR DATA ANALYSIS 

Several data sets of high precision are now publicly 
available. In addition to the COBE S data for small 
I there are data from BOOMERANG §, MAXIMA 
0, DASI (ll) and several other experiments {l|,[L3|. 
Wang, Tegmark and Zaldarriaga Jl2] (hereafter WTZ) 



have compiled a combined data set from all these avail- 
able data, including calibration errors. In order to test 
the robustness of our results, we do the analysis of N v 
for two different data sets. The first is the combined 
data of WTZ. The other consists of the COBE and 
Boomerang data, including the quoted calibration error 
of Boomerang ||. This second data set avoids possible 
systematics in the compiled data set. However, the final 
result for N v is practically the same for both data sets. 

The CMBR fluctuations are usually described in terms 
of the power spectrum, which is again expressed in terms 
of Ci coefficients as 1(1 + 1)C;, where 
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Ci = (\a lm \ 2 ). 
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The ai m coefBcients are given in terms of the actual tem- 
perature fluctuations as 



(4) 
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Given a set of experimental measurements, the likelihood 
function is 



£(6) ocexp ( -^[Ciey^x 
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where = fit,, Hq, n, r, . . .) is a vector describing the 
given point in parameter space, a; is a vector containing 
all the data points and C(8) is the data covariance ma- 
trix. This applies when the errors are Gaussian. If wc 
also assume that the errors are uncorrelated, this can be 
reduced to the simple expression, C oc e~ x / 2 , where 
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is a x 2 -statistics and iV max is the number of power spec- 
trum data points jl4| . In the present letter we use Eq. (||) 
for calculating \ 2 ■ 

The procedure is then to calculate the likelihood func- 
tion over the space of cosmological parameters. The ID 
likelihood function for N u is obtained by keeping N v fixed 
and maximizing C over the remaining parameter space. 

As free parameters in the likelihood analysis we use 
fi m , the matter density, Qb, the baryon density, Hq, the 
Hubble parameter, n, the scalar spectral index, r, the op- 
tical depth to reionization, and Q the overall normaliza- 
tion of the data. When large scale structure constraints 
are included we also use b, the normalization of the mat- 
ter power spectrum, as a free parameter. This means that 
we treat Q and b as free and uncorrelated parameters. 
This is very conservative and eliminates any possible sys- 
tematics involved in determining the bias parameter. We 
constrain the analysis to flat (f2 TO + fl\ = 1) models, and 
we assume that the tensor mode contribution is negligi- 
ble. These assumptions are compatible with analyses of 
the present data |12], and relaxing them do not have a 
big effect on the final results. For maximizing the likeli- 
hood function we use a simulated annealing method, as 
described in Ref. (lal. 



that 0.4 < h < 0.9 (h = H o /(100 kms" 1 Mpc" 1 )). 
For the ffo+BBN prior we use the constraint Hq = 
72 ± 8 kms" 1 Mpc -1 from the HST Hubble key project 
Jl6| (the constraint is added assuming a Gaussian distri- 
bution) and the constraint fl^h 2 = 0.020 ± 0.002 from 
BBN 0. Finally, in the # +BBN+LSS case, we add 
data from the PSC-z survey |1| to the data analysis. 

The neutrino density is to some extent degenerate with 
other parameters, particularly with the Hubble param- 
eter. Increasing the Hubble parameter allows for more 
neutrino species. In the same manner, decreasing n or 
ilbh 2 allows for more relativistic energy density. How- 
ever, N v is only slightly degenerate with these parame- 
ters. 

In Fig. 1 we show the likelihood functions for the two 
different data sets, assuming different priors. In the lower 
panels we show values of other parameters for the best 
fits. From this, it is evident that with only a weak prior 
on Hq, a large N v can be compensated by increasing Hq. 
As soon as the HST Hubble key project prior on H is 
added, the large values of N v are no longer allowed. 

From this figure it can also be seen that there is very 
little degeneracy between N v and n, il^/i 2 . Furthermore, 
the present data is entirely compatible with the BBN 
prior on f^/i 2 (as can also be seen in Fig. 1). Therefore, 
adding the BBN prior does not significantly change the 
analysis. 

In Table II, the best fit values and the 95% confidence 
limits on N v are shown for the two data sets, for different 
priors. Adding the prior h = 0.72 ± 0.08 from the HST 
key project gives a 2a upper limit of N v < 17 for the 
COBE+Boomerang data set and N v < 17.5 for the WTZ 
data set. 



B. LSS data 

Adding relativistic energy density also affects the mat- 
ter power spectrum because the growth factor on scales 
smaller than the horizon is decreased (see e.g. [|l9|). In 
Figs. 2 and 3 we show the CMBR and matter power spec- 
tra for the best fit models with different N v to the WTZ 
+ PSC-z data sets. 

Changing N v clearly also changes the matter spectrum, 
especially on scales of 0.01 — O.l/i Mpc -1 This fact can 
be used together with the CMBR data to constrain N v . 



A. Priors 



As was shown by Kneller et al. 0, different priors 
can significantly bias the derived confidence interval for 
N v . We therefore test the effect of different priors on 
the final result. Table I shows the different priors used. 
In the "weak" prior the only important constraint is 
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TABLE I. The different priors used in the analysis. 



prior type 




n b h 2 


h 


n 


T 


Q 


b 


"weak" 




0.008 - 0.040 


0.4-0.9 


0.66-1.34 


0-1 


free 




BBN + Ho 




0.020 ± 0.002 


0.72 ±0.08 


0.66-1.34 


0-1 


free 




BBN + H + LSS 




0.020 ± 0.002 


0.72 ± 0.08 


0.66-1.34 


0-1 


free 


free 
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FIG. 1. The top panels show the likelihood functions for the two different data sets, including different priors. The full lines 
are for the "weak" prior, the dotted for the Ho + fit/i 2 prior, and the dashed for the Ho + Qbh 2 +LSS prior. The lower panels 
show values of Ho, tlbh 2 and n for the best fit models. Horizontal full lines show the HST key project limit on Ho and the BBN 
prior on Qth 2 . 
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TABLE II. Best fit values and 2a (95%) limits on N v for 
different priors and the two different data sets. 



prior type 


WTZ 


COBE+Boomerang 


"weak" 


Q + H 


7 +l 7 > 


BBN + Ho 


8 +9 7 5 


4±f 


BBN + Ho + LSS 


6 +8 
°-4.5 


a -6.5 



6000 F 



Note that on even smaller scales, data from the Ly-a 
forest |20| can also be used. However, the very smallest 
scales are not so sensitive to N v because the shape of the 
power spectrum is not changed by adding radiation. The 
normalization is changed, but since we treat the overall 
normalization of the power spectrum as a free parame- 
ter, this will not have any effect. We therefore only use 
data from the PSC-z survey to give the LSS constraints. 
Adding the LSS data again tightens the constraint. The 
likelihood functions and best fit parameter values when 
LSS data is included can also be seen in Fig. 1. The 2a 
upper limits are now N v < 17 for the COBE+Boomerang 
data set and N v < 14 for the WTZ data set. 

For the WTZ data the upper bound is lowered from 
17.5 to 14 by adding LSS data. The effect can be seen 
in Figs. 2-3, for the model with N v = 14. Although this 
model can provide a very good fit to CMBR data, the 
shape of the matter spectrum becomes too shallow to 
obtain a decent fit. 

Very interestingly there is now also a non-trivial 
lower bound on N v which is N v > 2.5 for for the 
COBE+Boomerang data set and N„ > 1.5 for the WTZ 
data set. N v = is inconsistent with the data at roughly 
4er for COBE+Boomerang and 3a for WTZ. Indeed this 
result can be taken as the first real detection of the 
cosmological neutrino background at late epochs. From 
BBN considerations one already has the result N v > 2 
However, there is now an independent confirmation 
of the presence of relativistic energy density other than 
photons. Since the CMBR is only sensitive to radiation 
and not to the specific content, it is impossible to tell 
whether this radiation stems from the neutrinos as pre- 
dicted by the standard model, or from other light par- 
ticles. However, the standard result N v — 3.04 is in all 
cases compatible with the data at the 2a level. 

The incompatibility of N v — with data can also be 
seen in Figs. 2-3. Although a good fit to LSS data can 
be obtained, the fit to CMBR data is very poor. This 
is mainly because the first peak is too low due to the 
absence of the early integrated Sachs- Wolfe effect |2l] . 



III. DISCUSSION 

We have calculated bounds on the relativistic energy 
density present during recombination from the present 
CMBR and LSS data. The new data give a robust up- 
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FIG. 2. CMBR power spectra for the best fits to the 
WTZ+LSS data, for N u = (full line), 7 (dashed line) and 14 
(dotted line). The data points are from the WTZ compiled 
data set. 
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FIG. 3. matter power spectra for the best fits to the 
WTZ+LSS data, for N v = (full line), 7 (dashed line) and 
14 (dotted line). The normalization is arbitrary and the data 
points are from the PSC-z survey. 



per bound of N v < 17, but, perhaps more interesting, 
also give a lower bound of N v > 1.5/2.5 for the two 
different data sets analysed. Both bounds are interest- 
ing and non-trivial. Although the upper bound is much 
weaker than the bound N v < 4 found from BBN, it ap- 
plies to any type of relativistic energy density. The BBN 
bound can be avoided by putting some of the extra en- 
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ergj density in electron neutrinos, because these directly 



influence the neutron-proton conversion processes prior 
to BBN [|-|J2|]. The CMBR directly probes the energy 
density and is insensitive to flavour. The two constraints 
should therefore be seen as complementary. Furthermore, 
if there are massive particles decaying after BBN, but 
prior to recombination, the light decay products will add 
to the radiation density during recombination, but not 
during BBN. This is the case in some decaying neutrino 
scenarios ]T|,^3|], as well as in some scenarios with large 
extra dimensions pl| . 

The lower limit on N v is highly interesting because it 
provides the first strong indication of relativistic energy 
density other than photons around the epoch of recombi- 
nation. The value iV„ = is strongly disfavoured by the 
data, deviating from the best fit by 3<r for the WTZ+LSS 
data set and 4<r for the COBE+Boomerang+LSS data. 

Finally, although a central value higher than N v = 3 
seems to be preferred in the data, the standard model 
value N v = 3.04 is compatible with the present data at 
the 2cr level. This means that there is no significant indi- 
cation of non-standard physics contributing to N v at the 
recombination epoch. 
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